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We present results for form factors of semileptonic decays of D and B mesons in 2 + 1 flavor lattice QCD using 
the MILC gauge configurations. With an improved staggered action for light quarks, we successfully reduce the 
systematic error from the chiral extrapolation. The results for D decays are in agreement with experimental ones. 
The results for B decays are preliminary. Combining our results with experimental branching ratios, we then 
obtain the CKM matrix elements \Vcd\, \Vcs\, \Vcb\ and \Vub\- We also check CKM unitarity, for the first time, 
using only lattice QCD as the theoretical input. 



1. INTRODUCTION 

Semileptonic decays of B and D mesons play 
crucial roles in CKM phenomenology. The B de- 
cays such as B Tilv and B —f Dlv determine 
\Vub\ and \Vcb\, which are essential to constrain 
the CKM unitarity triangle. On the other hand, 
the D decays such as D nlv and D Klv pro- 
vide a good test of lattice calculations because 
corresponding CKM matrix elements \Vcd\ and 
\Vcs I are relatively well determined. In this paper, 
we report lattice calculations of semileptonic de- 
cays in unquenched (n/ — 2 + 1) QCD. By using a 
staggered- type fermion, which is fast to simulate. 
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Eigure 1. Result for CKM matrix. The first errors 
are theoretical, and the second experimental. 

for light quarks, we are able to reduce uncertain- 
ties from the "chiral" (m; mud) extrapolation. 
We calculate form factors for the above 4 different 
decays, from which the 4 CKM matrix elements 
are determined, as summarized in Eig. ^ The 
results for D decays are published in Ref. pp. 
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2. SIMULATION DETAILS 

We use Uf = 2 + 1 dynamical gauge configura- 
tions obtained with an improved staggered ( "Asq- 
tad") quark action on a lattice with « 1.6 
GeV, generated by the MILC collaboration 
For the valence light quarks we use the same 
staggered quark action, with the valence light 
quark [u, d) mass m™' equal to the dynamical 
light quark mass m^°^. The light quark masses 
we simulate range ^ <'mi< |ms, where uis is 
the strange quark mass. For the valence charm(c) 
and bottom(6) quarks we use a tadpole-improved 
clover action with the Fermilab interpretation . 
The hopping parameter for the c{b) quark is fixed 
from the Ds{Bs) mass. 

To form the heavy-light bilinears from the 
staggered-type light quark and the Wilson-type 
heavy quark, we convert the staggered-type quark 
to the naive- type quark, as in Refs. |4I5| . Rele- 
vant 3-point functions are then computed in the 
initial state meson rest frame using local sources 
and local sinks. We typically accumulate about 
500 configurations, and results at 2-4 source times 
are averaged to increase the statistics. 

For the matching factor of vector current Zy^ , 

we follow the method in Refs. jBu? , writing Zy — 
pv^iZy^ZyY^"^ . The flavor-conserving renormal- 

ization factors z'^^^^^ are determined nonpertur- 
batively from charge normalization conditions. 
For the remaining factor pv we use results in 
one-loop perturbation theory 

3. RESULTS 

3.1. D tt{K) and S -> tt 

The heavy-to-light decay amplitudes are pa- 
rameterized as 

{P\V^^\H) = f+{q^){pH +PP- AT + /o(g')A^ 

with q ~ ph — Pp, = (m|j — m^p)q'^/q^, 
V = Ph/itt-h, P± — Pp — Ev and E = Ep. The 
differential decay rate dT /dq^ is proportional to 
\ycKM\^\f+{q^)\^ ■ Below we briefly describe our 
analysis procedure; see Ref. ^ for details. 

We first extract the form factors f\\ and as 
in Ref. Q , and carry out the chiral extrapolation 
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Figure 2. mj-dependence and chiral fits for ff^'^. 

in mi for them at fixed E. To this end, we inter- 
polate and extrapolate the results for /|| and f± to 
common values of E using the parametrization of 
Becirevic and Kaidalov (BK) (Qj . We perform the 
chiral extrapolation using the NLO correction in 
staggered chiral perturbation theory (SxPT) 10\ 
We try various fit forms pP, as shown in Fig. [SJ 
and the differences between the fits are taken as 
associated systematic errors. 

We then convert the results for f± and /|| at 
mi = ruud, to /+ and /q. To extend /+ and /o 
to functions of q^ , we again make a fit using BK 
parameterization [H], 

f+i'f') ~ T^ ^2^77 /o(9"') ^ "i — ^~=TTq^ 

(1 — q'^)(l — aq-^) 1 — g^/p 

where q^ = q'^/rrij^,. We obtain 

ff"" = 0.23(2), a^'' = 0.63(5), P^"" = 1.18(5), 

for the B ^ TT decay, and 

f^"" = 0.64(3), a^"" = 0.44(4), /3^^ = 1.41(6), 

Z^^"^' = 0.73(3), a^^" = 0.50(4), /3^^ = 1.31(7), 

for the D decays, where the errors are statistical 
only. To estimate the error from BK parame- 
terization, we also make an alternative analysis, 
where we perform a 2-dimensional polynomial fit 
in [mi , E). A comparison between the two anal- 
yses are shown in Fig. |21 

Finally we determine the CKM matrix ele- 
ments (Fig. nj by integrating |/+((7^)p over 
and using experimental branching ratios jllll2 |. 
For \Vub\ we use the branching ratio for q"^ > 16 
GeV^ in Ref. ^^I- The systematic errors are sum- 
marized in Table ^ The results for D decays 
agree with experimental results 
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Figure 3. B ^ tt form factors from BK-based 
(filled) and non-BK-based (open) analyses. 



3.2. B- 

The B 



D 



D amplitude is parameterized as 

(DIVIB) = ^mBTUD X 

\h+{w){v + v'Y + h-(w)(v - v')% 

where v = vbI^ib-, v' = pd/tud and w = v ■ v' . 
The differential decay rate of i? ^ Dlv is pro- 
portional to the square of ^-'{w) , which is a linear 
combination of h+{w) and h-{w). We calculate 
the form factors at w = 1 by employing the dou- 
ble ratio method The light quark mass de- 
pendence for .^-"(1) is shown in Fig. 0] Extrapo- 
lating the result linearly to mi 0, we obtain 

7^+^(1) = 1.074(18)(16), (1) 



where the first error is statistical, and the second 
is systematic summarized in Tabled The system- 
atic error associated with finite lattice spacing is 
estimated by doing quenched calculations at dif- 
ferent lattice spacings and using different quark 
actions, and found to be small. 

Using Eq. and an experimental result for 
\Vcb\^{i) [111, we obtain \Vcb\ as given in Fig.[l| 

Table 1 

Systematic errors. 



decay 


D tt{K) 


B ^TT 


B D 


3-pt function 


3% 


3% 


1% 


BK fit 


2% 


4% 




mi extrap 


3%(2%) 


4% 


1% 


matching 


<1% 


1% 


1% 


a uncertainty 


1% 


1% 




finite a error 


9% 


9% 


<1% 


total 


10% 


11% 


2% 
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Figure 4. mj-dependence for J-b^d{^)- 

Since we have all 3 elements of the second row of 
CKM matrix, we are able to check a CKM uni- 
tarity using only our results as theoretical inputs; 

(iKdP + \Vcs? + \Vcb?f'^ - 1.00(10)(2). 
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